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Abstract 

We introduce quantized bipartite walks, compute their spectra, generalize the algorithms 
of Grover [G96| and Ambainis [Amb03] and interpret them as quantum walks with memory. 
We compare the performance of walk based classical and quantum algorithms and show that 
the latter run much quicker in general. Let P be a symmetric Markov chain with transition 
probabilities P[i,j], G [n]). Some elements of the state space are marked. We are 
promised that the set of marked elements has size either zero or at least en. The goal is to find 
out with great certainty which of the above two cases holds. Our model is a black box that 
can answer certain yes/no questions and can generate random elements picked from certain 
distributions. More specifically, by request the black box can give us a uniformly distributed 
random element for the cost of pQ. Also, when "inserting" an element i into the black box 
we can obtain a random element j, where j is distributed according to P[i,j]. The cost of 
the latter operation is p\. Finally, we can use the black box to test if an element i is marked, 
and this costs us p2. If S is the eigenvalue gap of P, there is a simple classical algorithm with 
cost O(po + (pi + p2)/Se) that solves the above promise problem. (The algorithm is efficient 
if po is much larger than p\ + p2.) In contrast, we show that for the "quantized" version of 
the algorithm it costs only O(po + (pi + P2)/ VSe) to solve the problem. We refer to this as 
the vfe rule. Among the technical contributions we give a formula for the spectrum of the 
product of two general reflections. 

1 Introduction 

The recent algorithm of Ambainis HAmb 03 1 which exploits properties of quantum walks on ver- 
sions of Johnson graphs opens up a new avenue for walk based algorithm designs. It claims as 
its direct predecessor the algorithm of Grover [G96| for database search. Ambainis's algorithm, 
like HG96I . is a search algorithm, but for many problems to which it applies, such as element dis- 
tinctness, it runs quicker than any application of HG96H . In the present article we give a thorough 
analysis of the type of walks that power the construction in [ Amb03 1, and use them in a very gen- 
eral algorithmic scheme. In our efforts to simplify the proofs we modify the paradigm that Grover 
has set up. 
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and in part by the National Security Agency (NS A) and Advanced Research and Development Activity (ARDA) under 
Army Research Office (ARO) contract number DAAD19-01-1-0506. 
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Ambainis's idea that a quantum walk can accelerate black box search for pairs of colliding 
elements demonstrates that diffusion-based quantum walks are more versatile than it had been 
thought previously. The consequences are numerous. Based on it Magniez, Santha and Szegedy 
[MSS) give a quicker algorithm for the oracle version of the triangle finding problem. Childs and 
Eisenberg [CE| streamline the analysis of HAmb 03l. and give further examples to its use. These 
papers, even the latter one do not change the proof structure of [ Amb03 1 significantly. 

In contrast, our discussion departs from that of HAmb0 3 1 at several points: 1 . We treat not only 
the Johnson graphs but all Markov chains; 2. We show how to eliminate the restriction Ambainis 
imposes on the state of the machine: we no longer need to request that it stays in a constant 
dimensional subspace of the entire state space. 3. We circumvent relevant parts of the proof by 
introducing "memory" for walks. As a consequence of I., 2., and 3. we obtain: 

Let P be a symmetric Markov chain with transition probabilities P[i, j], (i,j e [n]). Some 
elements of the state space are marked. We are promised that the set of marked elements has size 
either zero or at least en. The goal is to find out with great certainty which of the above two cases 
holds. Our model is a black box that can answer certain yes/no questions and can generate random 
elements picked from certain distributions. More specifically, by request the black box can give 
us a uniformly distributed random element for the cost of p . Also, when "inserting" an element 
i into the black box we can obtain a random element j, where j is distributed according to P[i, j}. 
The cost of the latter operation is p\. Finally, we can use the black box to test if an element i is 
marked, and this costs us If S is the eigenvalue gap of P, there is a simple classical algorithm 
with cost O(po + (pi + p2)/8e) that solves the above promise problem. (The algorithm is efficient 
if po is much larger than pi + p 2 .) In contrast, we show that for the "quantized" version of the 
algorithm it costs only O(po + (pi + p2) / V5e) to solve the problem. We refer to this as the \/~5e 
rule. Among the technical contributions we give a formula for the spectrum of the product of two 
general reflections. 

We settle at a walk model, directly derived from [Amb03|, that we call "bipartite." We focus 
only on walks that do diffusion for coin flip. Our formulas for the spectra and the eigenvalues of 
these walks should be useful in other contexts as well. 

In the quantum walk literature one can find two separate directions. Discrete time walks were 
introduced by Y. Aharonov, L. Davidovich, and N. Zagury lADZI and re-introduced by D. A. 
Meyer [ |Mey| . The properties of these walks were studied in one dimension by Ambainis, E. Bach, 
A. Nayak, A. Vishwanath, and J. Watrous: [ABNV|, and in general by D. Aharonov, A. Ambainis, 
J. Kempe, U. Vazirani IIAAKVI . 

Continuous time walks were introduced by E. Farhi, S. Gutmann HFGI . and they were studied 
by J. Roland, N. Cerf HRCI . by Wim van Dam, Michele Mosca, Umesh V. Vazirani [DMVl, by 
Childs and Goldstone ICGII and many others. Andrew M. Childs, Richard Cleve, Enrico Deotto, 
Edward Farhi, Sam Gutmann, Daniel A. Spielman HCCDFI show an example where a continuous 
quantum walk exponentially quicker traverses a graph than its deterministic counterpart. Continu- 
ous walks are intimately related to the paradigm of adiabatic computation. 

In this article we are concerned about special discrete time walks. The results are self- 
contained. 
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2 Notations 



We develop notations for bipartite objects, both classical and quantum. The symbol [n\ denotes 
a left-set of size n and symbol |to] denotes a right-set of size to. A classical bipartite walk takes 
place on the node set [n\ U |to]: 





Left 


Right 


Set: 


[n\ 


to] 


Vector: 


a = ...,a n \ 


6=|A,--.,/5m). 



We can concatenate left and right vectors into vectors of length n + to: 

(a,b) = (an, • • • ,a n ,/?i, . . ., A„). 

Left and right vectors and their concatenations are classical objects. To quantize them we form 
tensor products from them that lie in the Hilbert space C [ri|x|m] = C [n| <g> C |m] . 

Among the other not entirely standard notations we need are o, which means point-wise product 
of vectors and matrices and the square root of a vector or a matrix, which means point-wise square 
root: _ 

>/M[i,j] = y/M\iJ\. 

We apply the latter operation only if the elements of the matrix or vector are non-negative. Let 
Vi, . . . , Vk be vectors in the same Hilbert space. We denote the Gram matrix of v i, . . . , Vk (the 
matrix made of the inner products (v s , ft)) by Gram({vj}f =1 ). 

Linear operators will be denoted by any of lower case, upper case or Greek letters. The matrix 
of transition probabilities of a Markov chain will usually be denoted by P. The typical state space 
for our Markov chains is [n\ = {1, . . . , n}. Although we avoided using i as a running index in 
formulas containing the complex root of —1, the reader needs to use judgment about the meaning 
of % in each formula. 

3 Bipartite Walks 

Let [n\ and \m] be two disjoint sets of size n and m, respectively. A bipartite walk on [n\ U |to] 
is a probabilistic map with domain [n\ and range \m], and another probabilistic map with domain 
\m\ and range [n\. A probabilistic map is a stochastic matrix with rows and columns indexed by 
the elements of the domain and range of the map, respectively. Thus a walk is described by a 
pair of stochastic matrices (c, r) of dimensions (n, m) and (m, n), respectively. Since c and r are 
stochastic, we have that c[i,j] > and r[j, i] > for all 1 < i < n, 1 < j < m, Y^=i c [hj] = 1 
for all 1 < % < n, ^"=i r b> = ^ ^ or a ^ * — J ' — m - 

To implement the walk we need registers Left and Right, that hold values from [n\ and 
| to], respectively. An instantiation of the walk starts at an initial value a E [n\. We alternately 
set Right to j with probability c[Left, j] and then Left to % with probability r[i, Right], etc. If 
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y 

1/3 



i 



r C Prob = 1/2 xl x 1/2 x 1 x 1/2 

Walk Instantiation of length 5 

Figure 1 : Example to a bipartite walk, and one of its instantiations 

a comes from an initial distribution on [n\, then executing the walk for an even number of steps 
results in a distribution on [n\, and executing the walk for an odd number of steps results in a 
distribution on \m\. 

The quantized version of the walk takes place in the n x m dimensional Hilbert space of 
quantum registers |Lef t) | Right ) . The roles of c and r are taken up by diffusion operators 2C — I 
and 2R — I. Operator 2C — I is controlled by the first register and acts on the second register. 
Let q, (1 < i < n) be the vector of probabilities, where Cj[j] is the probability that c takes the i th 
element of [n\ into the j th element of \m\, and let rj, (1 < j < to) be the vector of probabilities, 
where rj[i] is the probability that r takes the j th element of \m] into the i th element of [n\. When 
|Lef t) is set to basis vector \i) we define the controlled diffusion operator on register |Right) as 

C 2\fCi\fCi — I\m}- 

The controlled diffusion operator of the |Lef t) register, when j is the value of the |Right) register 
is similarly defined by 

2 V ¥ jV ¥ f - hn\- 

Our focus of interest will be the operator, which does a controlled diffusion on | Right) and then 
another one on |Lef t). We redefine the above using the ket formalism. Let 

m 

Ci = ^2c[iJ]\i)\j), (l<i<n); 

n 

Writing all vectors in the basis \i) \j) (1 < % < n, 1 < j < to) we define the projection operators 

n 

c = J2V^V^ T ; (i) 



i=i 

m 



R = ( 2 ) 

3=1 
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Definition 1. The quantized version of the bipartite walk (c, r) is the pair (2C — I, 2R — I) of 
diffusion operators on C^ x ^. The two-step walk operator is p = (2R — I)(2C — I). 



In Section[T2]we shall compute the eigenvalues and eigenvectors of p in a more general setting. 
Equipped with these expressions in Section[8]we analyze the running time of our generalization of 
Grover/Ambainis type algorithms. 



4 The Discriminant Matrix 

The spectrum of p of the previous Section is completely determined by matrices c and r. But how? 
The following definition takes us one step closer to answering this question. 

Definition 2 (Discriminant Matrix). We define the discriminant matrix of a bipartite walk (c, r) 
as: 



M 





\/r o c T 



cor 



(3) 



The significance of the discriminant matrix will become clear from our Spectral Theorem. 
Below we give an equivalent definition of it, that will yield itself to a natural generalization. Notice 
that: 



Furthermore: 



(Vci | \/cr) = for all 1 < i < %' < n; 

(a/tj I y/ry) = 5jj/ for all 1 < j ' < f < m. 



(4) 
(5) 



Thus we can also express M as: 

M = Gram( A /cT, . . . , y/c^, i/F[, y/F^) - I n+m . 



5 A Spectral Theorem 

Let H be a Hilbert space and A be a subspace of H. The unitary operator that that leaves A 
invariant and takes all vectors in A 1 to their opposite is called a general reflection, and it is denoted 
by ref.4. Let A, B < H be defined via two separate orthogonal bases of unit vectors: 

A = (vi,...,v n ); 
B = (wi,...,w m ). 
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In Section[T2]we show the easy fact that for C = Ym=i v i v i* an d R = Y^jLi w j w j* tne operators 
2C — I and 2R — I are exactly ref^ and refg. Expressions ©, © and relations @ and © tie 
the problem of computing the spectrum and eigenvectors of p in Definition [T] to the problem of 
computing the spectrum of operator p = refgref^. 

Definition 3. We call A ± nB ± the idle subspace, and its orthogonal complement the busy subspace. 

The justification for the above definition is that the idle subspace lies in the kernel of operators 
C, R, and hence p acts in the idle subspace as the identity. Thus in order to get the spectral 
decomposition of p it is enough to compute it restricted on the busy subspace. 

Definition 4. For a = (a 1; . . . , a n | G (>l and b = . . . , (3 m ) e we define 

n 

a = ctiVi 
i=i 

rn 

i=i 

n m 

(a, 6)~ = a + b — y~] aiVj + } j PjWj. 

i=i j=i 

Every vector in the busy subspace has a convenient (albeit not unique) expression using the 
tilde: 

{A 1 n = (A,B)=A + B = {a + b\ ae C H ,b e C |m] }. 
We also need: 

Definition 5 (Discriminant Matrix (generalized)). The Discriminant Matrix of an ordered pair 
({vi, . . . , v n } , {wi, . . . , w m }) of orthonormal systems is: 

M = Gram(vi, . . . , v„, Wi, . . . , w m ) - / 

Our Spectral Theorem relates the spectrum and eigenvectors of p on the busy subspace to the 
spectrum and eigenvectors of M. The eigenvalues of M are symmetric to and are in the [—1, 1] 
range. 

Theorem 1 (Spectral Theorem). Let v 1; . . . , v n and w 1; . . . , w m be two orthogonal systems of 
unit vectors spanning spaces A and B, respectively. The eigenvectors and eigenvalues of the 
operator p = ref^ref^ on A + B are derived from those of the discriminant matrix M of the pair 
({vi}ie[n|, {Wj}je\m\) as follows: 
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Figure 2: The transformation that takes the eigenvalues of M into the eigenvalues of \x 
folds the [—1,1] interval and then expands it to the complex unit circle (dotted circle). 



Eigenvalue 


Eigenvector/Space 


1 


All vectors in An B. This space has dimension d\, where d\ is 
dimension of the eigen-space of M associated with eigenvalue 
1. 


2A 2 - 1 - 2iX\/l - A 2 


All vectors of the form a — \ b+ iy/l — A 2 6, where (a, b) is an 
eigenvector of M with eigenvalue 1 > A > 0. 


2A 2 - 1 + UXy/l - A 2 


All vectors of the form a — X b — i\J\ — X 2 b, where (a, b) is an 
eigenvector of M with eigenvalue 1 > A > 0. 


-1 


All vectors of the form a or b, where (a, b) is an eigenvector of 
M with eigenvalue 0. (This implies that (a, 0) and (0, b) are 
eigenvectors of M with eigenvalue too.) 



We prove the Spectral Theorem and other unproven claims of this section in Section [12| Let 
(a, b) G C' n ' u ' m l When {vi} ie [ n | and {wi} 3 - e | m ] specialize to {a/q}; and {y/fj}j, the tilde operator 
takes the following form: 



n I m \ m / n \ 

(a,&r = x>|i> X)M[i,j]|j> )+J2& E^»iw b*> = 
i=i \j=i j j=i \i=i j 

n m 
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The busy subspace is this case is the set of all vectors that can be expressed as above, and the 
discriminant matrix specializes to ©. 



6 Bipartite Walks from Ordinary Markov Chains 

A Markov chain with state set [n] is an n by n stochastic matrix P. We associate the (classical) 
bipartite walk (P, P) with chain P. The discriminant matrix of the walk is 

where D = \J P o P T . We call D the half discriminant matrix. Since M = ^ ^ ^ ^ ® D, 

its eigenvectors are of the form (a, a) and (a, —a), where a is an eigenvector of D. Matrix D in 
general is not stochastic, and it can be the zero matrix, for instance when P is the matrix of a cyclic 
permutation. A stationary distribution of a Markov chain P is a probability distribution zu on [n] 
such that zuP = zu. It is an eigenvector of P with all non-negative components. A Markov chain 
P is symmetric if P = P T . In this case D = P and its matrix is doubly stochastic. The uniform 
distribution on [n] is stationary for symmetric Markov chains. In the quantized case: 

Lemma 1. Let P be a symmetric Markov chain on [n] and let fi be the quantized walk operator 



associated with P. Let a = (-!=, . . . , -4=1, b — I -4=, . . . , -4=). Then u = a = b is a unit vector and 

jlU = U. 



Proof. We have 



M = ^Ev^ = E^Ev^MK)b) = (7) 



1=1 i=l v jr' = l 



1< j j'<n 



i=i v i=i v j=i 



(9) 



From Lemma H] it follows that u is an eigenvector of /i with eigenvalue 1. Direct calculation or 
Theorem |5] gives \u\ = 1. □ 



7 Walks with "Memory" 

We first study the classical setting. Let P be a symmetric Markov chain with state set [n]. Assume 
furthermore that elements of a subset G of [n] are marked, where G is either the empty set or 
\G\ > en. In this section we consider the task of designing an efficient algorithm that differentiates 
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in between the above two cases. At our disposal there are three subroutines, each with different 
associated costs: 



Name of the routine 


Description 


Cost 


PickUnif orm ( ) 


Picks a random i E [n] distributed according to zu 


Po 


ApplyChain (i) 


Applies the randomized map P on input i E [n] 


pi 


IsMarked (i) 


Returns if i E [n] \ G, returns 1 if i E G 


P2 



Under the above conditions what is the price of completing the task? The optimal algorithm 
depends on the ratios of p , pi, p 2 . If Po, Pu P2 h ave the same magnitude we just repeatedly 
pick random elements of [n] using PickUnif orm ( ) and test them with IsMarked (i) until 
we either find a marked element or declare that G is empty. We need to do this 0{l/e) times to 
achieve small constant error probability. However, if pi and p2 are much smaller than p , then the 
following algorithm will perform better: 

Algorithm FindMarked(fC) (Classical): 

Set i = PickUnif orm ( ) ; 
Do K times { 

If ( IsMarked (i) == ) i = ApplyChain ( i ) ; 

} 

Output IsMarked (i); 



If G = 0, the output of the above algorithm is always 0. On the other hand, if the eigenvalue gap 
of P is 5, \G\ > en and K is 1000/fe, it is easy to see that i will almost certainly "converge" to 
an element of G making the output 1 with high probability. The associated cost is: 

Classical cost of detecting large G = p + 1000(pi + p 2 )/(Se). (10) 

We shall "quantize" Algorithm FindMarked (K) and show that in the quantum case the 
right hand term of Expression (ITOb is replaced with 1000(pi + p 2 ) / \f5e. The merit of the result is 
that both e and 5 get square rooted, while what (almost) trivially follows from Grover is only that 
e gets square rooted. Our quantum machine will have registers: 



Control register: \b), where b E {0, 1}; 
Walk registers : | i ) \j ) , where i , j E [nj. 

We explain the role of the control register later. The register pair is used to perform a 
quantized bipartite walk with [n| = |ra] = [nj. Our new twist is that instead of quantizing P, we 
quantize: 

, def f if i G use map P on i; 

' ifiEG then map i into itself with probability 1 . 
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We might say that P' "remembers" if it ever sees a marked element. If G = then P' = P. For 
the general case, express: 



Pi P 2 
PI Ps 



Then the matrix form of P' is: 



with coordinate-division 





[n]\G 


G 


[n]\G 


Pi 


P2 


G 


P T 


P3 



p' 



Pi P2 

/ 



Above I is the identity matrix. The half discriminant matrix of the associated bipartite walk is: 

D = 



Pi 
/ 

Let us denote by v the quantum walk associated with P' and by fi the quantum walk associated 
with P. Define 



u 



l<i ,j<n 



P[i,j] 



n 



By Lemma[T]and © we have that nu = u. Therefore ji K u = u for every integer K . Our algorithm 
will utilize: 



1 . If G = then v = //, so 



|(u + u)/2| = 1. 



2. If \G\ > en and K > 1000/Vfe then 



u+v u 



< 3/4. 



The above formula explains the role of the control register. At the start we split the computation 
into two branches: in one branch we leave u untouched while in the other we apply v K on it. The 

1 1 
1 -1 



splitting operator is PI = (= H x ) = ^= ( * * ). The algorithm: 



Algorithm FindMarked(-ft') (Quantum): 

Put the pair of walk registers into the state u; 
Apply H on the Control register; 

/* This puts the system into state ^\^) u + ^l-O^' */ 

Do K times { 

If ( the Control register is |1) ) { 
apply v on the pair of walk registers; 

} 

/* Now the machine will be at state ^|0)u + ^j=\\)v K u */ 
Apply H^ 1 on the Control register; 
Measure the final state: 

If ( b = 1 or i G G ) output 1, else output 0. 
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Next we give a pair of subroutines that put the walk register \i) \j) into the state v\i) These 
routines implement the operators 2C' — I and 2R' — I, where (2C — /, 2R' — I) is the quantization 
of Markov chain P'. At our disposal, like in the classical analogue, we only have operators 2C — I 
and 2R — I, where (2(7 — I,2R — I) is the quantization of Markov chain P. While P does 
not not depend on the marked subset, P' does. Thus our subroutines will need to use subroutine 
IsMarked(i). They will also use an ancilla register \a) named Marked, where a e {0, 1}. This 
register stays |0) in between applications of the subroutines and is not to be confused with the 
Control register. 

Algorithm PerturbedWalk, diffusion on j: 

N)li)|0) ~^ where g(i) = IsMarked (i) ; 

If g(i) = apply (2(7 — /) on the pair of walk registers; 

/* This takes \i)\j)\g(i)) for g(i)=0 into the state 

\i) (£ /eW iVnJ] \f)) \g(i)) - \m\g(i)) v 

\i) \j) \a) — * \i) \j) \a © g(i)) , where g(i) = IsMarked(i); 

/* This step ^forgets'' g(i), since a=g(i) at this stage */ 
Similarly: 

Algorithm PerturbedWalk, diffusion on i: 

I*)b')b0'))/ where g(j) = IsMarked(j); 
If g(j) — apply (2R — I) on the pair of walk registers; 
/* This takes \i)\j)\g(j)) for g(j)=0 into the state 

(E^wav^MIO) \j)\gU)) - \m\g(j)) */ 

\i)\j)\ a ) ~> \i)\j)\ a ® 9(j)) r where g(j) = IsMarked(j); 

/* This step ^forgets'' g(j), since a=g ( j ) at this stage */ 

In the quantized setting we need to reinterpret costs p and pi. The cost of putting the machine 
into state u is p ; The cost of applying operators 2(7 — 7 and 2R — I on the pair of walk registers 
is pi. With this pricing the total cost of the operation FindMarked(fT) is 

O(p + K(pi + p 2 )). 

8 The y/5e rule 

Theorem 2. Let P be a symmetric Markov chain with state set [n], and let G be a subset of [n] 
marked via operator IsMarked (i) . Let 5 be the eigenvalue gap of P. Then for a randomly 
picked K e [1,1000/^]: 

1. IfG — FindMarked(K) outputs with probability 1; 

2- If\G\>en FindMarked(K) outputs 1 with probability at least 1/1000. 
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(The probabilities are both over K and the output of FindMarked(K ).) 



Proof. Recall that state of the machine during the execution of FindMarked(fQ evolves as: 

7f |0)tI+ 7! 



|0)0 -> |0}w -> -^|0}w+^|l)w -> (11) 



'" |0)m+ -> (12) 



ti + ^TI ^ TI-^tt (13) 



The proof of the last transition: 



-^|0)u+il)^u 




Q[n\xn] 



2 |0)(n + ^ M ) + -|l)( M -z/ K n). (14) 



The last arrow of (fTTT) represents a transition (actually the only one) that depends on the input, 
since v, the quantized version of the modified walk, P', depends on what G is. If G = then 
v = \i (see previous section), and since u is left invariant under u, the final state is \0}u, and the 
output is with probability 1 . This proves the first part of theorem. 
In order to prove the second part we show that if\G\ = en, then 

u + v K u 



< I (15) 

or the left walk register contains an element of G with probability 1/16. In fact, if e > 1/8 the 
latter is the case. To see this notice that the last operation does not effect the walk register, thus we 
get the same measurement for the walk register as if we measured it in the previous step. Taking a 
look now at Formula (|14j). and the fact that u us "uniform," we can easily see the claim. Therefore 
in the sequel we shall assume that e < 1/8. 

We use Theorem[l]to compute the spectral decomposition of v and in turn to compute the effect 
of v K on u. Recall that v is the quantization of the chain P', so in order to apply Theorem [T] we 
first need to determine the eigenvalues and eigenvectors of the half discriminant matrix 



D = x/P 7 ^ = J). 



For % E [n] let e, be the unit vector that takes 1 on i and elsewhere. Because of the block structure 
of D, for alH £ G the vector is an eigenvector of D. The remaining eigenvectors are those that 
are eigenvectors of Pi augmented with zeros on the coordinates corresponding to the elements of 
G. 

Lemma 2. The spectral radius of Pi is at most 1 — Se/2. 
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Proof. It is an easy fact that there exists a non-negative unit vector p' such that p'P\ = X'p' and A' 
is the spectral radius of Pi. Define p to be the augmentation of p' with zeros on the coordinates 
corresponding to the elements of G, and define a = Ym=i 7s e *- Since P is doubly stochastic we 
have aP = a. Consider the spectral decomposition of p in the basis formed by the eigenvectors of 
P: 

p = aa + ak a k, 

k 

where {a k } U {a} is a complete set of eigenvectors for P. Since P is symmetric, its eigenvectors 
are orthogonal and we have a 2 + J2 k a 2 . — 1. For a vector v G we denote by v G and Vq its G 
and [n] \ G components. |G| > en is equivalent to: 



(ag, ajr) = (oq, a) < 1 - e. 



Since p G = 0: 



a 



(p, a) = (p, a G + o,q) = (p, Oq) < \p\\oq\ < VT 



Hence, J2 k a\ = 1 — a 2 > e. Let X k be the eigenvalue of P associated with a k . Then pP 
aa + J2 k a k X k a k . Since the eigenvalue gap of P is at least 8, we have 

\pP\ 2 = a 2 + J2^ k X 2 k < 

k 

k 
k 

2J ot 2 k — 5 2^ o?k — 1 — $ £ - 



a 2 



a 2 



On the other hand 

\pP\ 2 >\ P D\ 2 = \p'P l \ 2 = X' 2 
which implies A' 2 < 1 — 5e, X' < 1 — 5e/2, as needed. □ 

From the above lemma one can anticipate the proof of the theorem: Any eigenvector of v that 
comes from the "quantization" of a shrinking eigenvector of P' has eigenvalue e ld , where 9 is 
separated away from by at least a constant factor times \f5e. This we deduce from the previous 
lemma and Theoremd Hence, if K is a random number in the range [1, 1000/\/fe], the phase 
shift, when applying v K on any of the above eigenvectors is some non-zero constant times n on 
expectation. Among the technical details we need to work out the most compelling (although not 
very hard) is that the projection of u on the space spanned by the above eigenvectors is large. 
Define 



ie[n]\G;je[n] 
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From our assumption that \G\/n < | we obtain: 



(u,u')= 



P[i,j] n-|GI 7 



> -. 



n n 

i£[n]\G; j'G[n] 



Since v is unitary so is v K , therefore 

- u K uf = \v K {u - u')\ 2 = \u- u'\ 2 = \u\ 2 + \u'\ 2 - (u, v!) - (vf, u) < ~. 



From this 



u + v K u 



I / I K / 1 2 i I /|2 i l I K"/ /\ 1 2 — -L i r . K / 1 2 . 

< - It + V M H M — M H U M — It < - \U + V U H . 

- 4 1 1 4' 1 4' v 71 - 4 1 1 8 

Set 

def i i K i \1 

amp K = | it + z/ m | . 

We are done if we show that 

Lemma 3. The probability that for a random K G [1, 1000/v^fe] the value o/amp K is at most 3 
is at least 1/6. 

Proof. We need to set up the stage to use Theorem[T]for the operator v. Since v is the quantization 
of P', the busy subspace will be A' + £>', where A' is generated by 

v-^Ys^PJJ] for i G [n] \ G; 

fori G G. 



and £>' is generated by 

w j = ^v / ^M K)|J) for j G [n] \ G; 

ie[n\ 

for ieG. 

The space spanned by the vectors (i G G) is invariant under z/ (from Lemma HI for instance) 
and it is orthogonal to v! = J2ie[n]\G "7s Vi- Define A" = (vi | 1 < i < n) and £>" = (wj | 1 < 
j < n). Let Z be the orthogonal complement of the subspace generated by \i)\i) (i G G). We 
have: 

u' G Z (16) 
z/|z = (ref S //ref^«)|x. (17) 



The discriminant matrix associated with the ({vi}j e [ n ]\G, { w j}je[n]\G) pair is 

Mi = 



Pi 
Pi 
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It follows from Lemma that all eigenvalues of M t are less than 1 — Se/2 in absolute value (the 
eigenvalues of Mi and those of Pi coincide up to a sign). Since the action of v on Z is the same 
as that of ref,B»ref_4» and v! E Z, we decompose u' according to orthogonal unit eigenvectors of 

refB»ref.4»: 

i 

u' = y^7fcZfc. 

fe=l 

Since u' E A", it lies in the busy subspace of operators ref^» and refg». We use Theorem [T] to 
claim that each z k on the busy subspace has eigenvalue 



e ie k dg 2ul-l + 2iu k ^l-ul 
where cu k is some eigenvalue of Mi. Consequentially \cu k \ < 1 — 5e/2. Hence 

|cos0 fc | = 2u 2 k -l < 2(l-fe/2)-l = 1-Se. 
Therefore when representing angles in [— n, n]: 

\e k \>V6e forl<k<£. (18) 

We have 

v K u> = J2 lk e^ K z k 

k=l 

l 

u ' + ^ u ' = £ 7fc (l + e***)z fc 

k=l 

i 

amp K = \v! + u K uf = ^| 7fc | 2 |l + e^| 2 . (19) 

k=l 

In Equation[T5]we used the fact that unitary operators have orthogonal eigenvector-systems. We 
show that for an individual k the expected value of |1 + e l9kK \ 2 for a random K E [0, lOOO/v^fe] 
is close to 2. Let us denote 1000/ y/5e by N. Then: 

1 - 

iV|l + e !W2 

AT ^ ' 



N 

K=l 



1 N 

N ^ 



K=l 

N N 
K=l K=l 

ie k (N+i) _ l ie k (N+i) _ i 

2 A I < 

N(e ie * - 1) N(e- i9 >° - 1) ~ 
2 2 

< 2.5, 



N(e i0 * - 1) JV(e- <9 * - 1) 
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when N is 1000/ v^fe. The last inequality comes from Inequality (TT81) Let 1 — p be the proba- 
bility that amp x > 3 for a random K. Then its expectation is lower bounded by 3 — 3p. Since 
2.5 Yli=i = 2-5 is an upper bound on the expectation, we have p > 1/6 as needed. □ 

□ 



9 Consequences 

As a first consequence we reprove the result of Ambainis [ Amb03 1 with a stronger implication: 

Theorem 3. Let X, Y be finite sets, f : X — ► Y be an oracle function and let TZ C Y x Y be a 
binary relation known to us. For H <0 X we define f(H) = {f(i) \ i G H}. Define 

p(f, a) = The probability that TZ PI (f(H) x f(H)) ^ for a random set H with size \X\ a . 

Then there is a quantum query machine with oracle f that can differentiate in between the cases 
when p(f, a) = andp(f, a) > e that runs in time 0{\X\ a + IOOOa/P^V^)- 

Proof. Let k = [|^| Q ] , n = ('*') , and P be the Markov chain on all k subsets of X with transi- 
tions: the probability that H goes to H' is zero if the the symmetric difference of H and H' is not 
two, and k Qx\-k) omerw i se - We quantize this chain, but with the caveat that with H we keep track 
of f(H). A set H is marked if TZ D (H X H) ^ 0. Since at all times we update f(H), it does not 
cost us queries to find out if H is marked. Hence p 2 = 0. The cost to perform 2C — I or 2R — I 
is 2. Hence pi = 2. To put the walk registers into the position 

u = !, E \H,f(H))\H',fm) 

V nk (\ x \ - k) {H ^\= 2 

costs 0(|X| Q ) queries. Hence p = \X\ a . The probability that an item is marked is p(f, a), which 
is either zero or at least e. Finally, the eigenvalues of the above Markov chain (derived from the 
Johnson graph) are well known to be — ~ k^x\-k)^^ ^ or J = 0, 1, . . . , r?,. Thus its eigenvalue gap 

k(\X\-k) ~k(\X\-k) 1 

The theorem now follows from Theorem |2| □ 

Corollary 1. There is a quantum query machine running in time 0(\X\ 2 ^) that differentiates in 
between the cases when f(i) ^ f(i') for alii ^ i' G X and when there are i ^ il G X such that 

Proof. Let TZ be the equality relation and set a = |X| 2 / 3 . □ 
The second consequence of Theorem |2] is immediate: 

Theorem 4. Let G be d-regular expander on {1, ... , n}. Then there is a version of Grover search 
that runs in \Jn steps and all transitions are done along the edges ofG. 

A special case, where we perform Grover search along the edges of the hypercube was studied 
by Julia Kempe |K|. 
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Figure 3: Two rectangular unit vectors with inner product 1/(2^). 

11 Operator Notation in Linear Algebra 

This part of the paper builds on ordinary linear algebra. Correspondingly, every vector (unless 
transposed) is a row vector. Operators (as opposed to the first part) are acting on the right, and in 
general we have to reverse the direction of all formulas of the first part involving vectors, stars and 
operators. For inner product we keep the angular notation, but we separate with comas rather than 
bars. In particular, if v is a vector and M and N are linear operators then vMN is a row vector 
that we obtain by applying M and N on v in this order. Also: 

(v, v) = vv* is a scalar; 

v*v is a dim v by dim v matrix. 

In order to represent elements of C' n ' <g> C' m l graphically we introduce Rectangular vectors. 
These are vectors with index set [n\ x \m], drawn in an array format and delimited by double bars 
in order to differentiate them from operators (see figure above). Regardless of their rectangular 
shape, they are row vectors in the sense that operators act on them on the right. If v and w are two 
rectangular vectors of the same dimensions, we can take their scalar product or we can create an 
operator by writing v*w that acts on C' n ' (g> C' m l. A nice thing about rectangular vectors is that 
we can conveniently express the tilde operator with them. Let M be annbym matrix and iV be 
an m by n matrix. We obtain the rectangular vector Mj by replacing all entries of M with zero, 
except those ones in the i th row. Similarly, we obtain Nj by replacing all entries of N T with zero, 
except those ones in the j th column. When N and M are given and fixed, and a = (ai, . . . , et n \, 
b = . . .,/3 m ), then 

n 

a = J2 aiMi ' (20) 

i=X 

m 

b = 5>N,-; (21) 

(a,6)~ = 5 + b. (22) 

We use the rectangular vector notation only in Section [TH where where we work out a specific 
example. 
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12 Spectra of Products of two General Reflections 



In this section we prove the claims of Section We need the definitions of that section with the 
modification that in this linear algebra inspired part we write operators on the right. Recall that 
subspaces A,B<H are defined via two separate orthogonal bases of unit vectors: 

A = (vi,...,v n ), 
B = (wi,...,w m ), 

and 

n 

i=l 
rn 

R = ^W* Wi . 

3=1 

Lemma 4. C is an orthogonal projection to A and R is an orthogonal projection to B. Also: 

rei A = 2C-I; (23) 
ref B = 2R-I. (24) 

Proof. We prove that C is an orthogonal projection to A and (|2"31 . The other claims are analogous. 

For V/ (1 < I < n) we have 

n n 

V;C = V;V*Vi = 2j(Vi, Vj)Vi = V;. (25) 
i=l i=l 

For any u which is orthogonal to all v^s we have: 

n n 



i=l i=l 



Equations d23t and (l2Tb prove that C is an orthogonal projection to ^4. Also from d2"5t and d2Tt : 

v^(2C - /) = 2vi - Vi = v ; and u(2C - I) = -u, and <E1> follows. □ 

Before computing the eigenvalues/vectors of fi we study the discriminant matrix of Definition 
13 M has a blocked structure corresponding to the subdivision of its rows and columns to vs and 
ws. The two diagonal blocks are and the two off-diagonal blocks are transposed conjugates of 
each other. 

Lemma 5. The spectral norm of M is at most 1. Furthermore, if {a, b) = (oti, . . . , a n , Pi, ... , (3 m ) 
is an eigenvector of M with eigenvalue 1 then 

n m 
i=l j=l 
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Proof. Since M is hermitian, all its eigenvalues are real. Let (a, b) = («i, . . . , a n , Pi, ... , (3 m ) be 
any unit. Then 

n m J n m \ 

(a, 6)M(a, b)* = 2j2J2 w j> = 2 ( a * Vi ' Yl } ' (27) 

i=i j=i \ i=i i=i / 

Let EILi N 2 = ££=i = 1 - ?• Since | ^ " =1 aiVi| = Jq and | E7=i^ w jl = V 1 ~ ?» 
the right hand side of (l2Vb is at most 2y/ q{l — q) < 1 with equality only if g = 1 — q = \. 
This implies the first part of the lemma. From the above it also follows that the right hand side of 
Equation dZTt is one iff Yl7=i a * v i = Ej=i ftj w y ^ 

We also show that the eigenvalues of M are distributed symmetrically to zero. The following 
observation will prove useful in many contexts: 

Lemma 6. If (a, b)M = (a', b') for some a, a' G CH b, b' G CH, then 

(a,0)M= (0,6'), (6, 0)M = (a', 0). 

Proof. Since the diagonal blocks of M are zero, (a, 0)M is of the form (0, b") and (0, 6)M is of 
the form (a", 0), which implies (a, b)M = (a", b"). But (a, b)M = (a', b'), which gives a" = a', 
b" = b', as needed. □ 

Lemma 7. 7f (a, 6), a G C^, 6 G C' m l w an eigenvector of M with eigenvalue A (a, —6) is an 
eigenvector of M with eigenvalue —A. 

Proof. From (a,6)M = (Aa, A6) LemmaElgives that (a,0)M = (0, A6) and (0, b)M = (Aa, 0). 
Then (a, -6)M = (— Aa, A6) = — A(a, —6). □ 

We denote the eigenvalues of M by 1 > Ai > A2 > • • • > A m+n and with g?a the dimension of 
the eigen-space associated with eigenvalue A. 

Let us now undertake the task of computing the spectrum and eigenvalues of 

fx = ref^ref B = (2C - I)(2R - I). 

Instead of \i it will be slightly more convenient to analyze the operator k = ~ (fJ>—I) = 2CR—C—R 
and its action on the busy subspace. Recall the definition of the tilde operation from Section 
From Lemma HJ 



dC = o; (28) 
bR = 6. (29) 

The proof of the theorem is powered by the following relations: 



aC = ((a,0)M)~; (30) 
= ((0,6)M)~. (31) 



21 



Proof. We prove only Equation (BUI) , the proof of (|3"T1) is analogous. Because of linearity it 
is enough to prove d3*0l for the basis vectors. Let 1 < i < n be arbitrary and let a = 
(0, . . . , 0, 1, 0, . . . , 0| be the unit vector, with a 1 in the i th position. Then a = Vj. We have: 

mm m 

aR = J~] VjW*Wj = ^(Vj, Wj)wj = y]M[vi,Wj]wj. 

j=i j=i j=i 

Above M[vj, Wj) means the entry of M indexed by the row associated with Vj and by the column 
associated with Wj. Thus aR = 6, where 

6 = |M[vi,w 1 ],...,M[v i ,w m ]). 

But this b is exactly (a, 0)M. (More precisely, (a, 0)M = (0, 6).) □ 

First we look at the action of k on A D 

Lemma 8. We characterize An B different ways: 

1. A vector is in Af~)B if and only if it can be written both as a and b. 

2. Every vector in Ad B is an eigenvector of ' k with eigenvalue 0. ( Or, equivalently, an eigen- 
vector of [i with eigenvalue 1.) 

3. a = b iff (a, b)M = (a, b), i.e. (a, b) is an eigenvector of M with eigenvalue 1. 

Proof. 1. is true by definition. For 2. assume that a = b. Then bC = aC = a and aR = bR = b. 
Hence 

2aCR -aC - aR = 2aR - a -b = 2b - a -b = 0, 

since a — b. Next we prove item 3. If a — b then bC = aC = a, so by d30b also (0, b)M = (a, 0). 
Similarly, aR = b, so by (l3TT> (a,0)M = (0,6). Therefore (a,b)M = (a, 6). Conversely, if 
(a, b) — (cti, . . . , a n , Pi, . . . , P m ) is an eigenvector of M with eigenvalue 1, then by Lemma |5] 
5 = 6. □ 

Next we shall create eigenvectors for k from the eigenvectors for M with eigenvalue less than 
1. Let (a, 6), a 6 C^, 6 G C' m ^ be an eigenvector of M with eigenvalue A ^ 0, i.e. 



{a, b)M = (\a, A6). 



Then Lemma|6lgives: 



(a,0)M = (0,A6); (32) 
(0,6)M = (Aa,0). (33) 



Combining these with (l30l and (I3TT) we get: 
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aR = Xb; (34) 

bC = Xa; (35) 

We would like to find a (3 such that v = a + (3b is an eigenvector of k. We use (l34l-(l35t to 
compute vk: 



VK =(a + (3b) (2CR -C-R) = (36) 

2aCR-aC -aR + 2(3bCR- (3bC - (3bR = (37) 

2A6 - a - Xb + 2(3X 2 ~b - (3Xa - (3b = (38) 

(-1 -/3X)a+ (X + 2(3X 2 - (3)b. (39) 

We conclude that as long as 

A + 2/3A 2 - (3 = P(-l- (3X), (40) 
v is an eigenvector of k with eigenvalue —1 — (3X. Let us express (3 from (l40l : 

(3 2 X + 2/3A 2 + A = 0. (41) 

Solving the equation gives (3 = — X ± \/X 2 — 1. Considering that ^ = 2k + 1 we can now write 
down the eigenvectors and eigenvalues of k and /i = 2k + I that we obtain from eigenvectors and 
eigenvalues of M. We summarize the formulas in the following two tables. The first table refers 
to the case of A G (0, 1): 





Eigenvector 


Eigenvalue 


M 


(a,b) 
(a, -b) 


A 

-A 


K 


a-Xb+iVl-X 2 b 


A 2 -l - iXVl-X 2 


a- Xb - iy/1- X 2 b 


A 2 - 1+ iXy/l- X 2 


(1 


a-Xb+iVl-X 2 b 


2A 2 -l-2iAVl-A 2 


a- Xb - i^l- X 2 b 


2A 2 - l + 2i\y/l- X 2 



The second table refers to the case when A = 0. In this case (a, b)M = (0, 0) implies that 
(a, 0)M = (0, 0) and (0, b)M = (0, 0). Therefore the zero subspace of M decomposes into the 
direct sum of the (possibly O-dimensional) subspaces {(a, 0) | (a, 0)M = (0,0)} and {(0,6) | 
(0, b)M = (0, 0)}. Since Equation (|4*TT) in this special case holds with every (3, we obtain: 





Eigenvector 


Eigenvalue 


M 


(a,0) 







(0,6) 





K 


a 


-1 




6 


-1 




a 


-1 




6 


-1 
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We are left to show that we have found all n + m — dim [A H B) orthogonal eigenvectors of 
the busy subspace. In the above tables we lined up the eigenvectors of k with the eigenvectors 
of M from which they originate to suggest a one-one correspondence. Observe that among the 
eigenvalues the correspondence is established by 

(M-side) ± A < — ► A 2 - 1 ± i AVl - A 2 (/c-side) 

The numbers seem to match, since M has n + m — d\ — d_i = n + m — 2d\ eigenvectors 
with eigenvalues in the range (—1,1). These correspond to the the same number of eigenvectors 
for k in the busy subspace with non-zero eigenvalues. In addition, the busy subspace contains 
d\ independent eigenvectors with eigenvalue zero. We get a total of n + m — d\ eigenvectors. 
By dim(^4 + B) = dim^4 + dim£> — dim (.4. n B) and and Lemma[8]the dimension of the busy 
subspace is also n + m — d\. We cannot walk away from the task, however, of showing that no 
dependencies occur among the eigenvectors we constructed. Since eigen-spaces associated with 
different eigenvalues are orthogonal, is sufficient to show that 

Lemma 9. Let |A| < 1. Then the dimensions of the eigen-spaces of n associated with eigenvalues 
A 2 — 1 ± i A\A — A 2 are d\ (each). Also, the dimension of the —1-eigen-space of k is do. 

Proof. Let A E (0, 1) and 

S x = {( a ,b) |aGC [n| , 6eC H , (a, b)M = X(a, b)}. 

Let r + and r_ be the operators from S\ to C' n ' x ' m l defined by 

r + : (a, b) — > d — Xb + iyl — A 2 b. 
t_ : (a, b) — > a — \b — iyl — A 2 b. 

We need to show that the images of r + and t_ have dimension d\ = dim S\. We prove that the 
kernel of both r + and t_ are trivial. We show this only for t + , since the proof for t_ goes in the 
same way. Let us assume, contrary to the lemma, that for some a E C' n ', b E m \ (a, b) E S\, 
(a, b) 7^ (0, 0) we have (a, b)r + = a — (A — iyl — A 2 ) b = 0. In fact it is enough to show that the 
assumption implies a = or b = 0, since both imply the other. By Lemma[8](a, (A — i\J\ — A 2 ) b) 
is an eigenvector of M with eigenvalue 1, which implies (a, 0)M = (0, (A — — A 2 ) b). Since 
(a, b) E S\, we also have (a, 0)M = (0, Xb). Hence 

Xb = (x - 2v / T 3 A 2 ) b, 

which, since |A| < 1 can happen only if b = 0, a contradiction. Note that the proof works for the 
A = case too. □ 

13 Norms and Inner Products 

In this section we show how to compute norms of vectors in A + B, and in particular we compute 
the norms of the eigenvectors we obtained in the previous section. We show that r + and r_ are 
scalar product preserving up to a constant scaling factor, and determine this constant. This gives 
an alternative proof to Lemma |9l 
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Lemma 10. For any a, a' G CH b, b' G CH it holds that 

(a, a') = (a, a'); (42) 

(b,b>) = (b, &'); (43) 

(a, b) = (a, 0)M(0, b)*. (44) 

Furthermore, if (a, b), (a', 6') G SV 

(a, a') = (6, 6') (45) 

(a, 9) = X(a, a) (46) 

Proof. Indeed, for a = (ct;|, a' = 6 = |$), 6' = 

(a, a') = ajVjV*a^ = OjC^ = (a a'); 

l<i<n l<i<n 

(1,9) = J2 h".,^ E W = ( bb 'h 

(5,6) = E E a * v i w^ = a i (v i ,w i )^= (a, 0)M(0, &)*. 

l<i<n l<i<m 

Consider now any (a, 6), (a', V) G S A . We have (a, 0)M(0, 6')* = (0, A6)(0, 6')* = A(6, &')■ 

But also, (a, 0)M(0, b')* = (a,Q)((Q,b')M*)* = (a,0)((0,b')M)* = (o,0)(Ao',0)* = A(a, a'). 

By the virtue of (14*41) the above shows not only (1431 but also (|4*6T ). □ 



Using Lemma [TO 



(5 - A6 + zv / T 3 A I 6j ~/ _ XV + zv 7 ! 3 ^ 2 b 1 ) = 

(a, a') + (b, 9) - X(b, a') - X(a, 9). (47) 



From Lemma flOl (b , a') = (a', b) = \(a', a) = A (a, a'). By introducing 7 = (a, a') we can write 
Expression d47l) as 7(2 — 2A 2 ). 
On the other hand 

((a, b), (a', b')) = (a, a') + (b, b') = 2 7 . 
We conclude that the scaling factor is 1 — 7 2 i.e. for every v, w G 5a: 

(r + (t;),r + H) = (l-A 2 )(t;, W ). 
Similarly we obtain that for every v, w G S\\ 

(t_(v) 1 t_(w)} = (1-\ 2 )(v,w). 

In particular: 

Theorem 5. The eigenvectors of fi in Theorem\I\have the norm 

Vl - A 2 , if|(a,&)| = landO < A < I; 
1, if a E An B and \ a\ = 1. 
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Figure 4: The Markov chain associated with Grover's algorithm, where items are marked with 
probability p 



In this section we give an example to the use of Theorem [T] for a Markov chain associated with 
Grover's algorithm. Here we present the "concise version" of the chain which we call the Grover 
Chain. The Grover chain has two states: marked and unmarked. The "full version," where differ- 
ent items correspond to different states has similar analysis. (We do not give a precise mathematical 
justification of the fact that clumping together all marked items and all unmarked items in the way 
we do gives formulas similar to those coming from the analysis of Algorithm FindMarked(fO for 
the chain P = ^E, where E is the all one matrix. Our example is interesting on its own right even 
without this connection.) 

Assume that the probability that an item is marked is p. The chain corresponds to the classical 
(non-quantum) algorithm, where at each step we move to a random item, but when we find a 
marked item we never move away from it. The transition of this chain takes an unmarked item to 
an unmarked item with probability 1 — p and to a marked item with probability p. On the other 
hand marked items alway go into marked items with probability 1. Figure [14| shows the Markov 
chain and its associated bipartite maps, c and r. The pair describing the walk is (c, r), where 



With eigenvectors a\ = (1,0) and a 2 = (0, 1), and eigenvalues 1 — p and 1 respectively. Let 



14 An Example 




The half discriminant matrix of the walk is 




b\ = (1, 0), b 2 = (0, 1). Then the eigenvectors of the discriminant matrix M = ( 

(a 1; bi) with eigenvalue 1 — p; (a 1; — &i) with eigenvalue p — 1; 

(ci2, 62) with eigenvalue 1; (a 2 , — b 2 ) with eigenvalue —1. 
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Let v be the quantized version of this bipartite chain. By Theorem[T]the busy subspace of v has 
eigenvectors: d\ — bi + p bi ± i^2p — p 2 bi and d 2 . The latter coincides with b 2 . Representing 
these as rectangular vectors: 



v 3 



y/1 - p (p + i^/2p - p 2 
y/p (-1 + p + iy f 2p 
( 



Vp 





V-2 



y/l-p (p-l^ 2p-p 2 ) y/p 

y/p (-1 +p - i\j2p -p 2 ) 







One can compute that \vi\ 2 = |f 2 | 2 = Ap — 2p 2 . This is consistent with Theorem |3 since the 
squared norm of (ai,&i) is 2, and Theorem |5] implies that the squared norms of v% and f 3 can 
be obtained from the norm of (ai,&i) by scaling it with the factor 1 — (1 — p) 2 . Theorem [l] 
gives that the eigenvalues associated with v\, v 2 and v 3 are 1 — Ap + 2p 2 — 2i(l — p)\j2p — p 2 , 
1 — Ap + 2p 2 + 2i(l — p) \j2p — p 2 and 1. What makes Graver's algorithm work is that 1 — Ap + 
2p 2 ± 2z(l - p)\j2p-p 2 = e ±ld , where 9 G &(y/p). Our version of Graver's algorithm (Section 
13) uses the initial state 

1 -p v 7 ^! 1 - p) 

This is the state we can easily produce (in the query model without cost, in the circuit model with 
small cost). Notice that u has a large component in the space spanned by v\ and v 2 . Define 



u 



u 



V2 



-vi + 



v i\ 



V2 



-v 2 



V-2 



Vp 



Then 



(u , v!) = a/1 -p. 
This again, should not surprise us because of Equation (fToT) . We have 



u'p k 



—le 



V2 



ie 



-i6K 



V2 



;v 2 , 



V2\ 



and our analysis of Algorithm FindMarked(iT) can go on like in the previous section. The goal of 
the present article is exactly to show how to shortcut much of the calculations we have made in 
this section. 
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